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Abstract 

A tracking detector is proposed for lepton flavor violation experi- 
ments {fj, — > e conversion, \i — > e + 7, \i — » 3e) consisting of identical 
chambers which can be reconfigured to meet the requirements for all 
three experiments. A pattern recognition and track reconstruction 
procedure based on the Kalman filter technique is presented for this 
detector. 

The pattern recognition proceeds in two stages. At the first stage 
only hit straw tube center coordinates, without drift time informa- 
tion, are used to reduce the background to a manageable level. At the 
second stage the drift time information is incorporated and a determin- 
istic annealing filter is applied to reach the final level of background 
suppression. The final track momentum reconstruction is provided by 
a combinatorial drop filter which is effective in hit-to-track assignment. 

The momentum resolution of the tracker in measuring monochro- 
matic leptons is found to be a p = 0.17 and 0.26 MeV for the \i — > e 
conversion and [i + — > e + + 7 processes, respectively. The tracker re- 
construction resolution for the total scalar lepton momentum is a p = 
0.33 MeV for the \x — > 3e process. The obtained tracker resolutions 
allow an increase in sensitivity to the branching ratios for these pro- 
cesses by a few orders of magnitude over current experimental limits. 



1 Introduction 

The goal of this work is to demonstrate a pattern recognition and track 
reconstruction procedure for a tracking detector which provides sufficient 
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resolution to be used in new lepton flavor violation experiments sensitive to 
branching ratios a few order of magnitude lower than current experimental 
limits. The observation of one of the three processes \i — > e conversion, 
\i —> e + 7 and /j — > 3e would provide the first direct evidence for lepton 
flavor violation in the charged lepton sector [U E] and require new physics, 
beyond the Standard Model. To reach such high sensitivity a substantial 
improvement in muon beam intensity [3] is required. This leads to the need 
for new designs of tracking detectors. A straw tube tracker proposed in this 
work consists of a set of identical chambers. A simple reconfiguration of the 
chambers allows meeting requirements for all three lepton flavor violation 
experiments. 

The signature of the \x — > e conversion (fi~ + iV — > e~ + N) process is 
clear: a single monochromatic electron in the final state with energy close 
to the muon mass m^. The fi + — > e + + 7 process has in the final state a 
monochromatic positron and a photon with an energy equal to half the muon 
mass. It is assumed that an external trigger is provided for the tracker. 
Photon reconstruction is not considered in this article because this study 
is devoted to pattern recognition and track reconstruction of electrons and 
positrons. The ii + — > e + + e + + e~ process has in the final state an electron 
and two positrons with energies from to m M /2 and an average energy equal 
to one third of the muon mass. 

In this work a pattern recognition and tracker reconstruction study was 
done for the // — > e conversion process, in the presence of background. A 
tracker resolution study for charged particles in ji + — > e + + 7 and ji + — > 
e + +e + +e~ processes was done without background. Our analysis is based on 
a full GEANT3 [4] simulation taking into account individual straw structure. 

2 Tracker description 

Muons stopping in a target produce electrons or positrons depending on 
the physical process. These charged particles move in a uniform magnetic 
field following helical trajectories. A tracking detector is used to measure 
the particle trajectory and thereby determine its momentum. The proposed 
tracker (Fig. [1]) consists of chambers with straw tubes placed transverse 
to a uniform magnetic field, which is along the tracker axis. Each straw 
tube measures the drift time corresponding to the radial distance at the 
closest approach of a charged particle from the sense wire. The tracker is 
operated in vacuum. The high muon beam intensity forces a tracker design 
without matter in a cylindrical central zone, to let the beam pass through 
the tracker. To get more redundancy for very rare signal events it is required 
that the measured trajectory should have two full turns in the tracker. This 
requirement sets limit on the minimal tracker length of about 300 cm. 

A tracker plane consists of two trapezoidal straw tube chambers (see Fig. 
[I]). The planes are distributed uniformly over the length of the tracker. 
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Figure 1: Schematic drawing of the tracker design. 

Rotating consecutive planes by an angle of 30° gives an effective "stereo" of 
crossed directions for 12 different views. Each chamber consists of one layer 
of straw tubes (60 straws) of 5 mm diameter. The straws are assumed to 
have wall thickness 25 //m and are constructed of kapton. Isobutane (C4H10) 
gas is assumed to fill the tubes. The tracker design allows changing the 
central zone radius, by moving the chambers closer to or further from the 
central axis, and changing the number of chambers along the tracker axis. 
This is needed to meet the different requirements for registration of electrons 
or positrons in the final state of each lepton flavor violation experiment. 

Two tracker chamber layouts and magnetic field configurations are stud- 
ied: one for the \x — > e conversion process and the second for \x — > e + 7 and 
ji — > e + e + e processes. The first tracker layout contains 108 planes and 
has the central zone radius of 38 cm. The second tracker layout contains 54 
planes and has the central zone radius of 10 cm. The tracker length for both 
layouts is the same, 300 cm. In the two layouts the tracker is immersed in a 
magnetic field of 1 T and 0.5 T, respectively. 

3 Pattern recognition 

Track reconstruction in lepton flavor violation experiments faces a significant 
amount of background hits in a tracker because of the high beam intensity 
needed to detect rare processes. A two stage procedure was developed to 
provide the pattern recognition in the tracker and to suppress background 
hits. At the first stage of pattern recognition only information on the centers 
of hit straws is used. The result of this stage is a significant suppression of 
background hits. Also, an approximate helix fitted to the straw hit centers 
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is found and used in the next stage of pattern recognition. To improve 
the suppression of background hits the second stage of pattern recognition 
uses drift time. A deterministic annealing filter DAF [5l [6] is applied at 
the second stage. The DAF effectively suppresses backgrounds by dealing 
simultaneously with multiple competing hits. 

3.1 Pattern recognition without drift time information 

A detailed pattern recognition and tracker reconstruction study was done for 
monochromatic electrons of 105 MeV from the // — ► e conversion process in 
the presence of background hits. In a uniform magnetic field the trajectory 
of a charged particle is a helix described by 5 parameters. The construction 
of the tracker allows significant simplification of the initial step of the pat- 
tern recognition by reducing the problem to two dimensional tracker views. 
Twelve tracker views are formed by planes through the tracker axis normal 
to the straws. In a typical event, a two-dimensional projection of the helical 
trajectory, a sine curve, is observed in a few views. There are approximately 
10 hits in each view, and 30 hits in the average in the event. The average 
number of background hits is 300 per event, which corresponds to an average 
straw hit rate of about 800 kHz. The hits are grouped in lobes (see Fig. 
EJ) with a typical gap between lobes of about 70 cm. As seen in Fig. [2] the 
sensitive area of the tracker measures a small part of a track, but combining 
hits from two lobes significantly improves the precision of the momentum 
reconstruction. 




Figure 2: Lobes in the tracker view. Points show straw hit centers. The 
sinusoidal projection of a helix is shown in red. The starting point of the 
track is not shown in the figure. 

The sinusoidal projection of a helix in each view is described by four 
parameters: x' , z' , Pl, Pt, 
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x[ = x' Q + clPt cos( 




(1) 



where Pl and Pt are the longitudinal and transversal momenta. The 
constant a equals 1/(2.998 B), where B is the magnetic field measured in 
Tesla. The coordinate z is common to all projections. Parameters x' and z' 
are defined for each given view such that the particle is created at the target 
with coordinate (x' Q + aP T , z' ). 

The reconstruction algorithm starts from a single tracker view. Four hits 
from this view are used to give a set of equations for helix parameters in 
the coordinate system related to that view. The system of four equations is 
solved to get the parameters Pl, Pt, z' ,x' Q [7]. All possible four hit combi- 
nations are considered. Only the combinations that survive a cut-off in the 
particle momentum (±20% of the expected value) are retained for subsequent 
analysis. 

Adding a fifth hit from a different view than the four hit combination 
allows defining a helical 3-D trajectory for the five-hit combination. 

The following approach is used to select signal hits. If we select a combi- 
nation with five signal hits the defined helical trajectory should correlate in 
space with other signal hits. For further analysis a five-hit combination is se- 
lected if the helical trajectory matches at least 15 additional hits in the road 
(±1 straw). A hit is rejected if it is not in any selected five-hit combinations. 
Due to the strong spatial correlations between the signal hits in comparison 
with the un-correlated background hits, the number of background hits is 
reduced drastically by applying this road requirement. A fit is applied to re- 
construct an average helical trajectory on the basis of the list of the selected 
tracker hits. This fit does not take into account multiple scattering. 

At the first stage of pattern recognition procedure the overall background 
rejection factor is about 130. The electron momentum from the fi — » e con- 
version process can be reconstructed in this stage with a standard deviation 
<y p = 0.45 MeV corresponding to relative momentum resolution cr p /p = 0.5%. 

Figure [3] shows the multiple view superposition of tracker hits for the 
sample event with 29 signal and 260 background hits before and after the 
pattern recognition procedure. There are no missed signal hits and two 
surviving background hits in this case. 

3.2 Kalman filter application to track reconstruction 

The Kalman filter (KF) [SJ E] addresses the general problem of trying to 
estimate at different points (1 < k < n) the state x k of a discrete process 
that is governed by the linear stochastic difference equation 



with a measurement mk = HkXk + £k- In the system of equations (J2]) 
stochastic processes such as multiple scattering and bremsstrahlung are taken 



x k = Fk_iXk_i + w k _i 



(2) 
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Figure 3: Plot of signal (blue dots) and background (red dots) tracker hits 
before (left) and after the pattern recognition procedure without drift time 
(right). 

into account by the process noise w^. The measurement noise is represented 
by £k- Qk and Vk are process noise and measurement noise covariances, 
respectively. In the absence of the last term Eq. (T5]) is the standard equation of 
motion with a propagator Fk-i (transport matrix). The F matrix propagates 
the state vector on one measurement plane to the state vector on the next 
plane combining position information with directional information. 

For a particle moving in a uniform magnetic field (see Fig. [T|) one has 
to choose the state vector parameters, define the initial state vector and 
calculate the transport matrix F, the projection matrix H, and the noise 
matrix Q. The state vector can be chosen in the form x k = (x, y, t x , t y , 1/pl) 
where x,y are the track coordinates in the tracker system, and t x = p x /pL, 
t y = Py/pL define the track direction. The projection matrix is given by 
H = (cos a, sin a, 0, 0, 0), where a is a tracker view angle. 

Due to multiple scattering the absolute value of electron momentum re- 
mains unaffected, while the direction is changed. This deflection can be 
described using two orthogonal scattering angles, which are also orthogonal 
to the particle momentum [10J. In terms of these variables the noise matrix 



is given by 
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For the variance of the multiple scattering angle the well-known expres- 
sion is used 



< e 2 >= (p /p) 2 [l + 0m8\n(t/X R )}t/X R 
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where po = 13.6 MeV, Xr is a radiation length and t is a distance traveled 
by the particle inside a scatterer. Energy losses are taken into account by 

pf = p— < dE/dx > t. 

There are three types of operations to be performed in the analysis of 
a track. Prediction x^" 1 is the estimation of the "future" state vector at 
position u k" using all the "past" measurements up to and including u k — 1". 

Filtering x£ is the estimation of the state vector at position based 
upon all "past" and "present" measurements up to and including "fc". At 
each step the filtered \ 2 is calculated. The total x 2 °f the track is given by 
the sum of the x\ contributions for each plane. 

Smoothing x£ is the estimation of the "past" state vector at position 
"/c" based on all "n" measurements taken up to the present time. The filter 
runs backward in time updating all filtered state vectors on the basis of 
information from all n planes. The mathematical equations describing these 
operations are given in [HE]. 

The prediction and filtering are applied consecutively to all points. When 
the last point is reached the smoothing is applied to all previous points. 
The result of smoothing is a reconstruction of the state vector x£ which 
defines particle coordinates and momentum in each plane. The KF approach 
described above works for a single point in each plane. Background hits create 
multiple competitive point in each plane. In this different approach 

based on the KF will be applied as described below. 

3.3 Pattern recognition with drift time 

The second stage in the pattern recognition procedure uses the hits selected 
and fitted helix obtained in the first stage and adds drift time information 
for the selected hits. The pattern recognition procedure can be improved 
by taking into account the measured drift time t™~ eas which is related to the 
radial distance r at the closest approach to the straw wire. The errors (a) 
in radius measurements are taken to be 0.2 mm. This radius r carries an 
ambiguity as to whether the track passed left or right of the wire. Left and 
right points are extracted from the intersections of a normal to the helix 
through the straw center and the circle of the radius r. 

A deterministic annealing filter (DAF) [5J E] is applied to suppress back- 
ground hits further. The DAF is a Kalman filter with re-weighted obser- 
vations. The propagation part of DAF is identical to the standard Kalman 
filter. In addition to background hits a left-right ambiguity for given hit 
creates multiple competing points at each KF step. At the DAF step all 
competing points are assigned to a single layer with weights. The filtered es- 
timate (measurement update) x£ at layer k is calculated as a weighted mean 
of the prediction xj^ 1 and the observations m^, i = 1, 2, ...n k . 
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By taking a weighted mean of the filtered states at every layer a prediction 
for the state vector x£* along with its covariance matrix C£* is obtained, 
using all hits except the ones at layer k. Initially all assignment probabilities 
for the hits in each layer are set to be equal but based on the estimated 
state vector x£* and its covariance matrix, the assignment probabilities of all 
competing hits are then recalculated in the following way: 

pi ~ v?(mL; H k x£* , V k + H k CrH£) 

where (p is a multivariate Gaussian probability density, V k is the variance 
of the observations. If the probability falls below a certain threshold, the 
hit is considered as background and is excluded from the list of the hits 
assigned to the track. At the initial step we cannot be sure of calculated 
probabilities due to insufficient information for the filter. This problem is 
overcome by adopting a simulated annealing iterative procedure [5j E] which 
allows avoiding a local minimum and finding the global one corresponding 
to the minimum chi-square for the track. The annealing schedule is chosen 
for the iteration n in the form V n = V(l + %) where the annealing factor 
/ > 1 and constant C 1. This insures that the initial variance is well 
above the nominal value V of the observation error but the final one tends 
to V. After each iteration the assigned probabilities exceeding the threshold 
are normalized to 1 and used again as weights in the next iteration, and so 
on. The iterations generally are stopped if the relative change in chi-square 
is less than a corresponding control parameter (typically of the order 0.01). 
Since we are dealing with a stochastic process, the best result can be reached 
repeating the DAF procedure for a few different annealing factors f (1.4 and 
2) and then choosing the result corresponding to the minimum chi-square. 

The application of the DAF procedure mainly allows resolving the left- 
right ambiguity. The DAF chooses wrong left-right points for 6% of the 
tracker hits if the radial distance r at the closest approach to the straw wire 
is greater then 0.25 mm. This corresponds to an average two tracker hits with 
the wrong left-right assignment per event. The number of background hits 
remaining after the second stage averages 0.38 hits per event in comparison 
with the primary 300 hits on average. The final background suppression 
factor is thus 300/0.38 ~ 800. Some of the signal hits are lost due to the 
selection (0.8 hits or 2.7% per event on average). The momentum resolution 
of the tracker after application of the DAF procedure is 0.25 MeV for the fi — > 
e conversion process. The pattern recognition procedure described above is 
thus effective in the rejection of background hits and also provides a good 
starting point for the track reconstruction. However, further improvement in 
the left-right assignment is necessary to achieve better resolution. For this 
we use the procedure described in the following section. 
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4 Momentum reconstruction for the \i — > e 
conversion process 

The momentum reconstruction is based on the hits selected by the two stage 
pattern recognition procedure described above. For straw drift tubes there is 
a set of left and right point for each straw hit due to the left-right ambiguity. 
An extreme possibility would be to apply the Kalman filter to all possible 
2 N combinations to reconstruct a track with N hits and to choose the combi- 
nation with the best \ 2 ■ This pure combinatorial approach may be the best 
but it is computationally unfeasible. The other extreme could be to apply 
the Kalman filter at each step to left and right points and to select the point 
providing the best x 2 . However this algorithm often selects wrong points. 

It has been proposed in the literature [5J E] to use a Gaussian sum filter 
(GSF) for solving the assignment problem in a track detector with ambi- 
guities [11]. In the GSF approach [12], both process noise and observation 
errors are modeled by a mixture of Gaussian densities and in the resulting 
algorithm several Kalman filters run in parallel. However the application 
of GSF as a momentum reconstruction procedure to this tracker does not 
provide sufficiently good results because of large gaps in the measured tra- 
jectory. In the tracker hits are distributed approximately uniformly inside 
lobes in the z-direction, but there are large gaps about 70 cm between lobes 
due to the motion of the electron outside the sensitive area of the tracker. A 
Kalman filter prediction step corresponding to a gap often leads to signifi- 
cant deviations from the true hit. As a result, in the GSF procedure the true 
component often is suppressed so significantly that it can not be recovered 
by subsequent Kalman filter steps. Therefore we need a different approach 
to treat the left-right ambiguity in the tracker. 

A combinatorial drop filter (CDF) has been developed to improve the mo- 
mentum reconstruction and efficiency for the tracker. The CDF algorithm 
uses forced minimization of the number of combinations when it reaches some 
maximum. So in the CDF approach the number of combinations always oscil- 
lates between N min and N max . That means that the Kalman filter runs with 
paths increasing in number with the combinatorial growth of combinations, 
but each time the number of combinations reaches N max , only the N min com- 
binations with the best x 2 are retained. The CDF approach is illustrated by 
the graph N min ^N max (2— >8) in Fig. HI Each vertex represents either a left 
or right point. Components surviving after the drop are indicated by small 
(magenta) circles. In the actual strategy the number of retained components 
was chosen to be 8 and the maximum number of components to be 32. 

The procedure begins with the first eight straw hits and builds 2 8 (256) 
possible hit combinations corresponding to left-right points. The Kalman 
filter forward and backward procedures are applied to these combinations. 
Only a few combinations which satisfy a rather loose \ 2 cu t, are retained. 
This step provides a certain precision for state-vector components. Each 
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Figure 4: Graph 2— >8 for CDF procedure. 

retained combination containing 8 points provides a starting point for the 
CDF. 

The Kalman filter reconstructs a trajectory of a particle in three dimen- 
sions. The trajectory is bent each time it crosses a tracker plane due to 
multiple scattering. Therefore, the reconstructed track is a set of helices 
that intersect at the planes. Fig. [5] shows the 3D trajectory reconstructed 
for a sample event. 



3D HELIX RECONSTRUCTED BY KF & 3D TRACKER HITS 




Figure 5: The 3D trajectory reconstructed for a sample event. 

On this scale the trajectory looks like a helix, but it consists of many 
helix parts. Fig. E] (left) shows the transverse projection of the trajectory 
for the sample event. In this projection the trajectory looks approximately 
like a circle. However if a specific region of Fig. [6] (left) is magnified one can 
see in Fig. [6] (right) that the shape of the circle is distorted due to multiple 
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3D HELIX RECONSTRUCTED BY KF S 3D TRACKER HITS 




scattering and energy loss. Three arcs of the trajectory are clearly seen in 
the Figure. 

There are two main tracker reconstruction features: momentum and angle 
resolutions. The momentum resolution is defined from the distribution of 
the difference between simulated and reconstructed momentum. The angle 
resolution is defined from the distribution of the angle between simulated 
and reconstructed momentum. 

In this study 10 5 events were simulated and reconstructed. The CDF 
procedure chooses wrong left-right points for 4% of the tracker hits (to be 
compared with 6% for the DAF) if the radial distance r at the closest ap- 
proach to the straw wire is greater then 0.25 mm. The distribution of the 
difference between the initial momentum reconstructed by the Kalman filter 
and the simulated initial momentum is shown in Fig. [7] (left) for the /i — > e 
conversion process. According to this distribution the intrinsic tracker reso- 
lution is a = 0.17 MeV if one fits the distribution by a Gaussian. Note that 
this tracker resolution is by about 30% better than the resolution provided 
by the DAF procedure. 

The distribution of the angle between the reconstructed and simulated 
momenta is shown in Fig. [7] (right). The most probable value of the angle 
distribution is 3 mrad. A calculation gives an average angle of multiple 
scattering in a straw to be 1.8 mrad for the fi — > e conversion process. This 
is consistent with the expectation that the angular resolution is defined by 
multiple scattering in a single tracker straw. 
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Figure 7: Reconstruction for fi — > e conversion process. The distribution 
of the difference between the reconstructed and simulated momentum (left). 
The distribution of the angle between the reconstructed and simulated mo- 
mentum (right). 



4.1 Momentum reconstruction for the fi + — > e + +7 pro- 
cess 

As mentioned in the introduction, this study and the one in the next section 
are less complete than that carried out for the \x — > e conversion process. No 
background hits were added to simulated events. The second tracker layout 
and field configuration were used for the momentum reconstruction study of 
the /i + — > e + + 7 process without background. The tracker is placed in a 
magnetic field of 0.5T which is half that used for the \x — > e conversion setup. 
The /i —>■ e conversion and the fi + — * e + + 7 processes have monochromatic 
charged leptons in the final state whose energies differ by a factor of 2. Due 
to the twofold change in the magnetic field the charged lepton trajectories for 
the fi + — > e + +7 and the fi — * e conversion processes are similar. The pattern 
recognition and momentum reconstruction procedures described above were 
applied to positrons of the fi + — > e + + 7 process. The distribution of the 
difference between the initial momentum reconstructed by the Kalman filter 
and the simulated initial momentum is shown in Fig. [8] (left). Fitting this 
distribution to a Gaussian gives an intrinsic tracker resolution of a = 0.26 
MeV. The distribution of the angle between the reconstructed and simulated 
momentum is shown in Fig. [8] (right). The most probable value of the 
angle distribution is 5 mrad. A calculation gives an average angle of multiple 
scattering in a straw to be 3.6 mrad for a positron in the fi + — > e + + 7 
process. 
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Figure 8: Reconstruction for the /i + — > e + + 7 process: The distribution 
of the difference between the reconstructed and simulated momentum (left). 
The distribution of the angle between the reconstructed and simulated mo- 
mentum (right). 



4.2 Momentum reconstruction for the /i — > 3e process 

The same setup as for the fi + — > e + + 7 process was used for a momentum 
reconstruction study of the fi + — ► e + +e + +e~ process. This process has in the 
final state an electron and two positrons with energies from to and an 
average energy equal to one third of the muon mass. The pattern recognition 
and momentum reconstruction procedure applied for \x —> e conversion and 
fi + — > e + + 7 processes was based on a single track search algorithm. The 
same procedure was applied separately for each lepton track of the /x + — > 
e + + e + + e~. The observable quantity of the process, total scalar momentum 
of charged leptons p tot = J^Pi, was reconstructed. 

The distribution of the difference between reconstructed total scalar mo- 
mentum and muon mass is shown in Fig. [9] (left). Fitting this distribution 
by a Gaussian gives a tracker resolution for total scalar momentum of a = 
0.33 MeV. The distribution of the angle between the reconstructed lepton 
momentum and the simulated momentum is shown in Fig. [9] (right). The 
most probable value of the angle distribution is 7 mrad. A calculation gives 
an average angle of multiple scattering in a straw to be 5.4 mrad for charged 
lepton momentum equals m^/3. 



5 Conclusion 

A pattern recognition and momentum reconstruction procedure based on a 
Kalman filter technique for a straw tube tracker was developed to suppress 
background and reconstruct track momentum for the lepton flavor violating 
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Figure 9: Reconstruction for the /i — > 3e process: The distribution of the 
difference between reconstructed total scalar momentum and muon mass 
(left). The distribution of the angle between the reconstructed and simulated 
momentum (right). 

processes \i — > e conversion, fi + —>■ e + + 7 and fi — > 3e. The simple modular 
construction of the tracker allows meeting requirements for lepton momentum 
measurements for all lepton flavor violation processes mentioned above. 

The pattern recognition procedure allows suppressing the background hits 
by a factor of 800. Only 3% of signal hits are lost by this procedure. 

The momentum resolution of the tracker to register monochromatic lep- 
tons was found to be a p =0.17 and 0.26 MeV for the /i — > e conversion with 
background and the /z + — > e + + 7 processes without background, respec- 
tively. The tracker resolution for the total scalar lepton momentum is a p = 
0.33 MeV for the fi — > 3e process without background. For the /1 — > e con- 
version, /i + — > e + + 7 and /i — *■ 3e processes the most probable values of an 
angle between the simulated and reconstructed momenta (3, 5 and 7 mrad) 
are in a good agreement with multiple scattering angles in a single tracker 
straw. The obtained tracker resolutions allow an increase in sensitivity to 
the branching ratios for these processes by a few orders of magnitude over 
current experimental limits. 

We wish to thank A. Mincer and P.Nemethy for fruitful discussions and 
helpful remarks. This work has been supported by the National Science 
Foundation under grants PHY 0428662, PHY 0514425 and PHY 0629419. 
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